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ABSTRACT 

In this paper we provide a short overview of the scope and strong future potential of a multi-messenger 
approach to gravitational-wave astronomy, that seeks to optimally combine gravtitational wave and elec¬ 
tromagnetic observations. We highlight the importance of a multi-messenger approach for detecting grav¬ 
itational wave sources, and also describe some ways in which joint gravitational wave and electromagnetic 
observations can improve the estimation of source parameters and better inform our understanding of the 
sources themselves - thus enhancing their potential as probes of astrophysics and cosmology. 

Subject headings: gravitational waves, multi-messenger astronomy, electromagnetic follow-ups 


1. Introduction 

As we mark the centenary of General Relativity, the 
rapidly emerging field of gravitational wave astronomy 
stands on the threshold of a new era, with the global 
network of so-called ‘second generation’ ground- 
based interferometers preparing to begin operations 
(The Virgo Collaboration 2009; Harry & LIGO Scientific 
2010). The approximately ten-fold improvement in 
sensitivity of these detectors compared with the first 
generation LIGO and Virgo detectors (Aasi & et al. 
2015) means that the first direct detections of gravi¬ 
tational waves (GWs) are eagerly anticipated within 
the next few years - with the most likely sources 
expected to be the inspiral and merger of com¬ 
pact binary systems: neutron star-neutron star (NS- 
NS) binaries; black hole-black hole (BH-BH) bi¬ 
naries or NS-BH binaries. The challenges asso¬ 
ciated with detecting such sources, and exploit¬ 
ing them as probes of astrophysics and cosmology, 
have attracted considerable effort in recent years - 
including thorough investigations of the GW de¬ 
tection rates that can be expected (Abadie & et al. 
2010), the most likely electromagnetic (EM) coun¬ 
terparts that will accompany them (Metzger & Berger 
2012) and the efficiency with which those EM coun¬ 
terparts could be detected (Nissanke et al. 2013b; 
Cowperthwaite & Berger 2015). What is clear from 


^Xilong.Fan@glasgow.ac.uk 


this work is that significant benefits will derive from 
joint observations of both GW signals and their EM 
counterparts - giving rise to a new observational field 
that has been termed ‘multi-messenger’ astronomy. 

The benefits of a multi-messenger approach are 
manifest for both the detection and characterisation 
of GW sources. Perhaps the most immediate benefit 
° ^is°a'precise identification of the source sky location 
- alleviating the relatively poor sky localisation that 
is possible from ground-based GW observations alone 
(Wen et al. 2008; Wen & Chen 2010; Nissanke et al. 
2011; Schutz 2011; Eairhurst 2011; Sidery et al. 2014). 
Another clear benefit derives from the potential iden¬ 
tification of an EM ‘trigger’, or associated emission, 
that may streamline and enhance GW searches by re¬ 
ducing the number of free parameters, or at least sig¬ 
nificantly restrict the range of their values that must 
be explored (Mohanty et al. 2004; Abbott et al. 2008). 
Equally, a multi-messenger approach may involve us¬ 
ing GW ‘triggers’ to prompt targetted EM searches 
for the counterpart of the GW source (Abadie & et al. 
2012a; Aasi & et al. 2014) - a successful outcome of 
which would indeed both identify a precise sky loca¬ 
tion and permit follow-up observation of the source’s 
host galaxy, allowing for example measurement of its 
redshift. More generally, joint GW-EM observations 
may in turn lead to a better and more robust infer¬ 
ence of the source parameters, thus opening the way 
to using GW sources more effectively as astrophysical 
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probes. 

There are very significant observational and com¬ 
putational challenges associated with these multi¬ 
messenger approaches, however. From a GW point 
of view, there is the challenge of carrying out a suf¬ 
ficiently rapid real-time analysis of the interferometer 
data - thus allowing meaningful information to be sup¬ 
plied to the EM community quickly enough to permit 
follow-up observations of any prompt EM emission as¬ 
sociated with the GW event (Abadie & et al. 2012a,b; 
Evans et al. 2012; Aasi & et al. 2014). To act upon the 
results of such an analysis then presents major chal¬ 
lenges from an EM point of view: e.g. understanding 
what are the EM counterparts themselves, and what 
are their signatures across the EM spectrum, and then 
searching efficiently for those signatures - the latter 
task rendered difficult in view of the poor sky localisa¬ 
tion provided by the GW triggers, particularly during 
the first few years of operation of the second genera¬ 
tion network (Singer et al. 2014). 

In this paper we will consider some of these multi¬ 
messenger issues in more detail for the particular case 
of NS-NS or NS-BH compact binary coalescences 
(CBCs). We will very briefly summarise our cur¬ 
rent understanding of the expected EM counterparts 
of these sources and highlight some of the observa¬ 
tional issues this raises for their detection through 
rapid follow-up observations across the EM spectrum. 
We will also briefly discuss how GW searches can pro¬ 
vide suitable triggers to the EM community on an ap¬ 
propriately short timescale. Einally, we consider the 
astrophysical potential for multi-messenger observa¬ 
tions of CBCs - highlighting one specific cosmolog¬ 
ical question that may be addressed in the fairly near 
future with these data. 

Of course a multi-messenger approach is poten¬ 
tially crucial for all GW sources - not just CBCs - 
and may become increasingly important as the ‘de¬ 
tection era’ of gravitational-wave astronomy unfolds. 
Eor example, looking ahead to possible third gen¬ 
eration ground-based detectors such as the proposed 
European Einstein Telescope (Punturo & et al. 2010) 
there are excellent prospects for multi-messenger sci¬ 
ence (involving not just GW and EM observations 
but also cosmic rays and neutrinos, see for exam¬ 
ple (Aartsen & et al. 2014)) with a wide range of as¬ 
trophysical sources. The interested reader can find 
more details in e.g. Chassande-Mottin et al. (2011); 
Andersson et al. (2013) but we do not consider the 
broader context for future multi-messenger astronomy 


any further in this paper. 

2. Joint detection of GW and EM emission: con¬ 
straints from the event timescale 

To make best use of a multi-messenger approach to 
observing GW sources requires both a solid theoret¬ 
ical understanding of the signatures produced by the 
EM counterpart of the source, together with the abil¬ 
ity to observe these signatures straightforwardly. Con¬ 
sider the example of a NS-NS or NS-BH CBC merger 
event, for which the likely counterpart is thought to be 
a short-duration gamma ray burst (sGRB) (Paczynski 
1986; Berger etal. 2005). Metzger & Berger (2012) 
have addressed in detail the question of what might be 
the most promising electromagnetic counterpart of an 
sGRB, and they identify four key features (which they 
term ‘cardinal virtues’) that such a counterpart should 
possess. It should: 

1. be detectable with present or upcoming tele¬ 
scope facilities, provided a reasonable allocation 
of resources 

2. accompany a high fraction of GW events 

3. be unambiguously identifiable, such that it can 
be distinguished from other astrophysical tran¬ 
sients 

4. allow for a determination of ~arcsecond sky po¬ 
sition 

We will return to virtues 2, 3 and 4 in the next sec¬ 
tion, but the physical nature of a sGRB event makes 
satisfying virtue 1 already a challenge - particularly 
when one is considering EM searches prompted by 
a GW trigger - due to the very short timescales in¬ 
volved. According to the standard picture, in the first 
one or two seconds immediately after the merger the 
sGRB results from rapid accretion onto a centrifugally 
supported disk that surrounds the merged object; this 
accretion powers a highly relativistic jet that, due to 
beaming, is only seen as a sGRB within a narrow half¬ 
opening angle of a few degrees (Burrows et al. 2006; 
Soderberget al. 2006). Interaction of the jet with the 
interstellar medium produces prompt afterglow emis¬ 
sion that is visible in X-rays on a timescale of sec¬ 
onds to minutes after the merger, and in the optical on 
a timescale of hours to days - but again only within 
a narrow solid angle around the jet. On the other 
hand it is believed that the event will also emit fainter 
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isotropic radiation, in the optical and infra-red, on a 
timescale of days after the merger; this emission is 
caused by the decay of heavy elements synthesised 
by the merger ejecta as it interacts with the surround¬ 
ing environment, and is referred to as a ‘kilonova’ 
(Li & Paczyhski 1998; Kulkarni 2005; Metzger et al. 
2010; Kasen et al. 2015). We will see in the next sec¬ 
tion that kilonovae may be key to the successful iden¬ 
tification of an EM counterpart for many CBC events. 

Thus, as is summarised in e.g. Singer et al. (2014) 
and Singer & Price (2015), the nature of the EM coun¬ 
terparts of a NS-NS or NS-BH merger, and their rapid 
evolution in time, requires that the GW trigger be 
generated as quickly as possible - and ideally on a 
timescale of seconds to (at most) hours - particularly 
if the event afterglow is to be observed in X-rays. 

Such a low latency, template-based data analysis 
pipeline has been developed for searching for CBC 
signals with the ground-based interferometer network 
- and the need for a very rapid trigger generation is 
one of several key factors driving this search. (Eor 
more details see Cannon et al. (2012) ). In the sixth 
LIGO science run and third Virgo science run, the 
inspiral search pipeline Multi-Band Template Anal¬ 
ysis (MBTA) (Beauville & et al. 2008) was able to 
achieve trigger-generation latencies of 2-5 minutes 
(See e.g. (Abadie & et al. 2012a)). Interesting co¬ 
incident triggers identified by the MBTA were sub¬ 
mitted to the Gravitational-wave Candidate Event 
Database (GraCEDb) ^ However, with the advanced 
detector network, given that sources will spend up to 
10 times longer within the detection band, the low- 
latency trigger generation of inspiral signals will be¬ 
come even more challenging. 

The slow evolution in frequency of an inspiral sig¬ 
nal allows one to reduce the sampling rate of both 
the data and templates. The highly similarity be¬ 
tween neighbouring templates in the filter banks al¬ 
lows one effectively to reduce the number of filters or 
to use so-called ‘parallel infinite impulse response’ fil¬ 
ters. We refer the reader to the LLOID (Low Latency 
Online Inspiral Detection) (Cannon et al. 2012) and 
SPUR (Summed Parallel Infinite Impulse Response) 
(Luan et al. 2012; Hooper et al. 2012; Liu et al. 2012) 
pipelines for details of the implementation of these ef¬ 
forts. 
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3. Joint detection of GW and EM emission: iden¬ 
tifying the source location 

Prom analysis of the GW interferometer network 
data alone, an approximate sky position for a CBC 
source can be obtained by exploiting the difference in 
arrival times of the GW signals at the different detector 
locations (Aasi & et al. 2013). In the first few years of 
operation of the advanced detectors the typical ‘error 
box’ derived in this manner will, unfortunately, have 
an area of several hundred square degrees (Singer et al. 
2014). The situation will improve substantially in later 
years as more interferometers are added to the global 
network, and by the early 2020s - with the anticipated 
inclusion of the proposed LIGO India detector - al¬ 
most 50% of binary NS merger systems detected by 
the network are expected to be localised to within 20 
square degrees on the sky. Looking even further ahead, 
recent work (Raffai et al. 2013; Hu et al. 2015) investi¬ 
gating the optimal sites of future, third generation, de¬ 
tectors has applied various metrics to define and quan¬ 
tify the performance of a given network configuration 
- including the ability of the network to localise the po¬ 
sition of a typical CBC source on the sky. However, a 
fundamental limitation is set by the scale of the Earth, 
and the maximum difference which this imposes on 
the time of arrival of GWs at well separated detectors. 
Thus, for any ground-based detector network, present 
or planned, we can expect that GW data alone will lo¬ 
cate many sources to a precision of only tens or even 
hundreds of square degrees. Clearly, then, identifying 
a unique EM counterpart would be a crucial benefit of 
a multi-messenger approach as it can provide a much 
more precise sky location for the GW source. 

So how might such an EM counterpart be identi¬ 
fied, given such relatively poor sky localisation infor¬ 
mation, in the context of the ‘cardinal virtues’ high¬ 
lighted above? As discussed in Metzger & Berger 
(2012), sGRBs are not in themselves likely candidates 
since the beaming of their relativistic jets will mean 
that most merger events will not be seen as a sGRB 
at the Earth - even if a sGRB does indeed ‘accom¬ 
pany a high fraction of events’. Similarly the prompt 
X-ray and optical afterglow emission associated with 
the sGRB may also not be seen at the Earth - al¬ 
though if such an afterglow were observed through 
rapid follow-up observations with a high-energy burst 
monitor satellite such as Swift, then this can in princi¬ 
ple determine a very precise (~ 10 arcseconds) sky lo¬ 
cation for the GW source (Evans et al. 2012). As con- 
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eluded in Metzger & Berger (2012), therefore, kilono- 
vae are perhaps the most promising candidates for the 
EM counterpart of a NS-NS or NS-BH merger event 
since their isotropic emission gives them the poten¬ 
tial to satisfy all four cardinal virtues. On the other 
hand their comparative faintness and their limited du¬ 
ration still suggests that their detection will be very 
challenging, involving rapid, wide-held EM follow-up 
campaigns. In fact to date only two candidate kilo- 
novae events have been observed (Tanviret al. 2013; 
Yang etal. 2015). 

Whatever the nature of its EM counterpart, a key 
step towards its study is the (hopefully unique) iden- 
tiheation of the host galaxy of the GW source. To 
this end we can develop a general framework for com¬ 
bining GW and EM information under the straight¬ 
forward assumption that there are common parame¬ 
ters - such as sky location and distance - among the 
observations made of a GW source, its EM counter¬ 
part and their host galaxy. Such a multi-messenger ap¬ 
proach has already proven fruitful even in the absence 
of any observed GWs. Eor example, searches for GWs 
associated with the gamma-ray bursts GRB070201 
(Abbott & et al. 2008) and GRB051103 (Abadie et al. 
2012) have ruled out the possibility that their progen¬ 
itors were binary neutron stars (NS) sources in An¬ 
dromeda galaxy (M31) and M81, respectively. 

One can consider the question, then, of how the de¬ 
tection efficiency of an EM counterpart might be im¬ 
proved by conducting wide-held EM follow-ups that 
make use of pre-existing galaxy catalogs. This issue 
has been investigated in Hanna et al. (2014); by tak¬ 
ing into account the GW measurement error in both 
distance and sky location they estimated that an aver¬ 
age of ~ 500 galaxies are located in a typical GW sky 
location error box for NS-NS mergers with Advanced 
LIGO (~20 deg^), up to range of 200 Mpc. The au¬ 
thors then found that the use of a complete reference 
galaxy catalog could improve the probability of suc¬ 
cessful identiheation of the host galaxy by ~ 10-300% 
(depending on the telescope held of view) relative to 
follow-up strategies that do not utilize such a catalog. 

Of course a more complete treatment ideally should 
involve a comprehensive end-to-end simulation that 
takes fully into account the characteristics of the global 
GW network, the spatial distribution and event rate 
of CBC events, the multi-wavelength light curves of 
their EM counterparts (together with the event rate and 
lightcurves of a range of possible ‘false positives’ with 
which these counterparts might be confused, bearing 


in mind virtue 3 above) and the sensitivies and observ¬ 
ing strategies adopted by the multi-wavelenth EM tele¬ 
scopes used to search for them. Much work remains to 
be done in this area, both on more detailed theoreti¬ 
cal modelling of the EM counterparts themselves (e.g. 
kilonovae) and on investigating the impact and efficacy 
of different follow-up observing strategies. Some ex¬ 
cellent progress has already been made, however, par¬ 
ticularly in Nissanke et al. (2013b) where the authors 
point out that: 

• dedicated Im class optical and near infra-red 
telescopes with very wide-held cameras are 
well-positioned for the challenge of hnding EM 
counterparts 

• a comprehensive catalog, out to z ~ 0.1, of fore¬ 
ground stellar sources, background active galac¬ 
tic nuclei and potential host galaxies in the lo¬ 
cal universe - and probing deeply (i.e. approx¬ 
imately 2 magnitudes deeper than the expected 
magnitudes of the EM counterparts) - is needed 
to assist identiheation of the EM counterparts 
amidst the larger numbers of false positives. 

To date no such all sky galaxy catalog is avail¬ 
able. Among the many existing galaxy catalogs, 
the Gravitational Wave Galaxy Catalog (White et al. 
201 l)(GWGC;(White et al. 2011) has been specih- 
cally compiled for current follow-up searches of opti¬ 
cal counterparts from GW triggers - although the au¬ 
thors acknowledged that the catalog is increasingly in¬ 
complete at larger distances. Nuttall & Sutton (2010) 
proposed a ranking statistic to identify the most likely 
GW host galaxy (drawn from the GWGC) based on 
galaxy distance and luminosity and the sky position 
error box. This ranking method has been adopted 
in the design of an EM follow-up pipeline ( such 
as Nuttall et al. (2013)) and follow-up observations 
(such as Aasi & et al. (2014); Kanneretal. (2012); 
Ando et al. (2013)). 

A novel Bayesian approach to identifying the GW 
host galaxy, incorporating astrophysical information 
has been proposed in Ean et al. (2014). A merit of this 
approach is that both the rank and the posterior proba¬ 
bility of a galaxy hosting the GW source are estimated 
with the help of GW-EM relation models. Using a sim¬ 
ulated population of CBCs the authors found that (i) 
about 8%, 4%, and 3% of injections had 50%, 90%, 
and 99% respectively of the probability to be included 
in the top 10 ranked galaxies in the GWGC, and (ii) the 
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first ranked galaxy had a 50% probability of being the 
true GW host galaxy in about 4% of injections. These 
results are dominated by the GWGC distance cut of 
100 Mpc, compared with the expected reach of ~200 
Mpc for Advanced LIGO and Advanced Virgo at de¬ 
sign sensitivity. The method could easily be extended 
and applied to deeper galaxy catalogs, however. 

A comprehensive, complete all sky catalog of all 
objects (galaxies, stellar sources and so on) may not be 
essential for a single GW-EM joint observation, since 
the GW sky map will only cover a small fraction of the 
whole sky. A “galaxy survey on the fly” approach has 
been proposed to reduce the required effort and thus 
enhance the immediate availability of useful catalogs 
(Bartos et al. 2015); in this approach a rapid galaxy 
survey using 1-2 m class telescopes can efficiently cat¬ 
alog those galaxies covered by one GW detection vol¬ 
ume within a short period of time. This rapidly com¬ 
piled galaxy catalog could then very quickly be pro¬ 
vided to other telescopes, to aid further electromag¬ 
netic follow-up observations of e.g. kilonovae, as well 
as other sources. 

Whatever detailed strategies are employed in the 
future, it is clear that to identify the common GW 
and EM source one needs multi-wavelength observa¬ 
tions capable of rapidly and efficiently targeting mul¬ 
tiple objects. A community of ‘multi-messenger’ as¬ 
tronomers - including the LIGO and Virgo collabo¬ 
rations and a significant and growing number of tele¬ 
scopes across the world - is now being assembled for 
exactly this purpose^^. 

4. Beyond detection 

The promise of multi-messenger astronomy - and 
in particular of joint GW and EM observations of 
CBCs and their EM counterparts - lies not just in en¬ 
hancing the prospects for detecting these events but 
also in improving the estimation of their parameters, 
and ultimately their use as astrophysical and cosmo¬ 
logical probes. 

Consider again, for example, the case of a CBC 
event. In the absence of any spin, the inspiral phase 
of a CBC in a circular orbit has a gravitational wave¬ 
form that depends on the following nine parameters: 
four intrinsic parameters - the two masses and two in¬ 
trinsic constants of integration, which define the phase 
evolution of the waveform; five extrinsic parameters 
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- luminosity distance, sky location, inclination angle 
(the angle between the orbital angular momentum and 
the line of sight) and GW polarization. All nine pa¬ 
rameters govern the amplitude of the signal that each 
GW detector should ‘see’. So we see immediately 
that the identification of a unique EM counterpart will 
begin to reduce the dimensionality and complexity of 
this parameter space, providing tight constraints on the 
sky location and the luminosity distance (e.g. via a 
measured redshift for the host galaxy - although de¬ 
riving a luminosity distance directly from this will be 
complicated by the effects of galaxy peculiar veloci¬ 
ties). Knowledge of the luminosity distance from an 
EM counterpart will also help to break parameter deg- 
neracies that may exist for the GW source, such as that 
between distance and inclination angle (Nissanke et al. 
2010). 

More generally, joint GW-EM observations should 
allow for better characterization of the signal progen¬ 
itor, and a richer interpretation of the results of GW 
searches - offering, for example, greater insight at the 
population level into the relationship between the GW 
source, its EM counterpart(s) and their progenitor en¬ 
vironment. An excellent example of such an approach 
is the relationship between CBCs and sGRBs that we 
have discussed in detail in this paper, but as the detec¬ 
tion era unfolds no doubt other similar examples will 
emerge. The insights derived from these joint analy¬ 
ses may improve our understanding of the evolution¬ 
ary processes that led to the GW event, and thus per¬ 
haps provide better prior information about the likely 
sites (e.g. in terms of host galaxy morphology, colour, 
metallicity etc) of future events. Moreover, as our un¬ 
derstanding of these deep, underlying astrophysical re¬ 
lationships improves, new data analysis methods will 
be built to explore them - combining not just GW 
and EM observations but also data from neutrino and 
cosmic-ray telescopes (Adrian-Martmez & et al. 2013; 
Aartsen & et al. 2014). 

Underpinning this joint approach at all stages is the 
fundamental assumption that these multi-messenger 
phenomena all emanate from the same object, and 
obey physical relationships that share some com¬ 
mon parameters. Bayesian inference then provides 
a natural framework in which to incorporate multi¬ 
messenger astrophysical information and optimally 
estimate those parameters. The method introduced in 
Ean et al. (2014, 2015) explores some specific exam¬ 
ples of this framework in action; these include reduc¬ 
ing the area of sky over which astronomical telescopes 
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must search for an EM counterpart, improving the in¬ 
ference on the inclination angle of a NS-NS binary and 
estimating the luminosity of a sGRB. 

What specific astrophysical questions will be tack¬ 
led via a multi-messenger approach in the future? We 
end by briefly considering one such example that is 
anticipated for the advanced detector era. 

It was recognised nearly thirty years ago by Schutz 
(1986) Schutz (1986) that CBCs have the potential 
to be precise luminosity distance indicators via mea¬ 
surement of the time-evolution of their amplitude, fre¬ 
quency and frequency derivative during the inspiral 
and merger phase. This has given rise to the idea of 
so-called gravitational-wave standard sirens, by anal¬ 
ogy with the standard candles of EM astronomy - al¬ 
though the CBCs do not require any assumption to be 
made about their intrinsic ‘luminosity’ (which essen¬ 
tially depends on the masses of the binary stars) as this 
can be inferred directly from the observed data at the 
same time as the luminosity distance. In more recent 
years there has, therefore, been much interest in the po¬ 
tential use of standard sirens as cosmological probes, 
via calibration of the luminosity distance redshift rela¬ 
tion. 

The expected reach of advanced detectors will 
be too shallow to permit exploration of dark en¬ 
ergy models and the accelerated expansion of the 
Universe - although such models could certainly be 
investigated by third generation detectors such as 
the Einstein Telescope (Sathyaprakash et al. 2010; 
Zhaoetal. 2011) or spaceborne missions such as 
eLISA (Holz & Hughes 2005; Deffayet & Menou 
2007; Amaro-Seoane & et al. 2012). However, a 
realistic target for the upcoming global network 
of advanced detectors is measurement of the Hub¬ 
ble constant, Hq, using standard sirens. Recently 
Nissanke et al. (2013a) have investigated the efficacy 
of such a measurement, and conclude that a preci¬ 
sion of about 1% on //o is possible from observations 
of about 30 NS-NS mergers within a few hundred 
Mpc - using the anticipated future global interefer- 
ometer network that includes KAGRA and LIGO In¬ 
dia. Such a precise value would certainly be com¬ 
petitive with the EM results expected on a similar 
timescale from e.g. the James Webb Space Telescope 
(Freedman & Madore 2010), and in any case would be 
an extremely useful adjoint to a purely EM determina¬ 
tion of Hq using the traditional cosmic distance ladder 
since it would be subject to a completely different set 
of systematic uncertainties. 


A standard siren measurement of Hq will present a 
major multi-messenger challenge for several reasons. 
Firstly, to estimate the Hubble constant of course re¬ 
quires comparison of distance with redshift, and the 
latter will not generally be measurable from GW data 
alone (but see also below for discussion of some inter¬ 
esting alternative approaches). Indeed the luminosity 
distance estimates for the sirens will in any case be 
degenerate with redshift because of the mass-redshift 
degeneracy in post-Newtonian CBC waveforms. By 
measuring the redshift of the siren’s host galaxy the 
degeneracy can immediately be broken. However, this 
measurement of course first requires the prompt obser¬ 
vation of an EM counterpart and the unique identifica¬ 
tion of the host galaxy - steps which will be subject 
to all of the multi-messenger issues discussed in the 
previous section. 

Notwithstanding these potential difficulties, and 
their resulting impact on the final error budget for Hq 
(for which the the estimate of 1% in Nissanke et al. 
(2013a) may therefore be somewhat too optmistic), 
the prospect of a gravitationally-calibrated value of 
the Hubble constant is nevertheless extremely excit¬ 
ing - and is likely to be one of the main targets for 
gravitational-wave astronomy over the next decade. 

It is interesting to note that some other approaches 
to breaking the mass-redshift degeneracy, and/or de¬ 
termining redshifts (and hence cosmological parame¬ 
ters) from GW data alone, have been proposed. For 
example Taylor et al. (2012) assume that there exists 
a universal (rest frame) mass distribution for NSs at 
different redshifts and by comparing the measured 
(redshifted) mass distribution of NSs with the local 
mass distribution show that one could infer statisti¬ 
cally the redshifts of the sources and hence derive in¬ 
directly the value of Hq. Their results show that, in 
this way, second generation interferometers should be 
able to infer the Hubble constant with ~ 10% accuracy 
from about 100 events. Their analysis is extended in 
Taylor & Gair (2012) to consider the cosmological po¬ 
tential of GW-only observations with third generation 
ground-based detectors. 

In a similar manner, various authors MacLeod & Hogan 
(2008), Petiteau et al. (2011) have proposed that the 
identification of the host galaxy - and thus the determi¬ 
nation of its redshift - may be carried out statistically, 
using prior information about the spatial distribution 
of the galaxies in the sky localisation error box pro¬ 
vided by the GW data alone. Using this formalism 
del Pozzo (2014) Del Pozzo (2014) suggests that the 
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Hubble constant could be determined to a few percent 
from observations of about 50 sirens with the advanced 
detector network. 

A third, highly promising, possibility for constrain¬ 
ing cosmological models using GW observations alone 
has been proposed by Messenger & Read (2012) 
Messenger & Read (2012), exploiting the effect of 
tidal deformations on NS-NS binary systems during 
the inspiral that provide additional contributions to the 
phase evolution of their gravitational waveforms. Re¬ 
cently their approach has been studied further, and 
is predicted to be capable of determining redshifts 

to a precision of 10-20% for GW sources in the 

local Universe observed with the Einstein Telescope 
Messenger et al. (2014). The potential of this approach 
for determining cosmological parameters with Ein¬ 
stein Telescope observations has been further explored 
in Del Pozzo et al. (2015). 

While these various methods suggests the intrigu¬ 
ing possibility, therefore, of measuring Hq - and in¬ 
deed other cosmological parameters - without using 
any EM observations directly, the generally lower pre¬ 
cision of these estimate does, nonetheless, underline 
that a multi-messenger approach will usually be more 
effective than a GW-only analysis. 

5. Summary 

In this paper we have discussed some of the ad¬ 
vantages - both for the detection of GW sources and 
also for the estimation of their parameters and their 
astrophysical exploitation - of a multi-messenger ap¬ 
proach that seeks to combine optimally GW and EM 
observations. Eocussing mainly on the inspiral and 
merger of NS-NS and NS-BH binaries, which are be¬ 
lieved to be the progenitors of short duration GRBs, 
we have also highlighted some of the important ob¬ 
servational challenges that need to be overcome in or¬ 
der that a multi-messenger approach may be fully ex¬ 
ploited. These challenges present significant logistical 
and computational constraints for the analysis of data 
from the ground-based network of advanced GW inter¬ 
ferometers. This is because of the very short timescale 
(from seconds to hours to days at most) associated 
with the EM counterparts of these events and the rel¬ 
atively poor sky localisation provided by the GW data 
alone - which will significantly complicate the search 
for a unique EM counterpart. Nevertheless substan¬ 
tial progress has already been made towards establish¬ 
ing a community of ‘multi-messenger’ astronomers. 


working closely together towards the goal of mak¬ 
ing joint GW-EM observations, and there are excel¬ 
lent prospects for this emerging new field over the next 
decade. 

Looking further ahead, the potential of multi¬ 
messenger astronomy would appear to grow even 
stronger with the possible advent of third genera¬ 
tion ground-based interferometers such as the Einstein 
Telescope and the possible launch of a spaceborne 
GW detector such as eLISA. Although we did not dis¬ 
cuss these future missions and projects in this paper, 
they should offer exciting science possibilities such as 
probing the equation of state and internal structure of 
neutron stars, constraining models of core-collapse su¬ 
pernovae and mapping the detailed structure of space- 
time around black holes. As general relativity enters 
its second century, the prospects for testing Einstein’s 
theory under extreme cosmic conditions using multi¬ 
messenger data look very bright. 
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Theoretical Physics in Beijing. 

REFERENCES 

Aartsen, M. G., & et al. 2014, Phys. Rev. D, 90, 
102002 

Aasi, J., & et al. 2013, ArXiv e-prints, 1304.0670 
—. 2014, ApJS, 211,7 

—. 2015, Classical and Quantum Gravity, 32, 074001 

Abadie, J., & et al. 2010, Classical and Quantum Grav¬ 
ity, 27, 173001 

—. 2012a, A&A, 541, A155 

—. 2012b, A&A, 539, A124 

Abadie, J., et al. 2012, ApJ, 755, 2 

Abbott, B., & et al. 2008, ApJ, 681, 1419 

Abbott, B., et al. 2008, ApJ, 683, L45 

Adrian-Martmez, S., & et al. 2013, J. Cosmology As- 
tropart. Phys., 6, 8 

Amaro-Seoane, & et al. 2012, Classical and Quantum 
Gravity, 29, 124016 


7 


Andersson, N., et al. 2013, Classical and Quantum 
Gravity, 30, 193002 

Ando, S., et al. 2013, Rev. Mod. Phys., 85, 1401 

Bartos, I., Crotts, A. P. S., & Marka, S. 2015, ApJ, 801, 
LI 

Beauville, R, & et al. 2008, Classical and Quantum 
Gravity, 25, 045001 

Berger, E., et al. 2005, Nature, 438, 988 

Burrows, D. N., et al. 2006, ApJ, 653, 468 

Cannon, K., et al. 2012, ApJ, 748, 136 

Chassande-Mottin, E., Hendry, M., Sutton, P. J., &. 
Marka, S. 2011, General Relativity and Gravitation, 
43, 437 

Cowperthwaite, P. S., & Berger, E. 2015, ArXiv e- 
prints, 1503.07869 

Deffayet, C., & Menou, K. 2007, ApJ, 668, L143 

Del Pozzo, W. 2014, Journal of Physics Conference 
Series, 484,012030 

Del Pozzo, W., Li, T. G. E, & Messenger, C. 2015, 
ArXiv e-prints, 1506.06590 

Evans, P. A., et al. 2012, ApJS, 203, 28 

Eairhurst, S. 2011, Classical and Quantum Gravity, 28, 
105021 

Pan, X., Messenger, C., & Heng, I. S. 2014, ApJ, 795, 
43 

—. 2015, Astrophysics and Space Science Proceed¬ 
ings, 40, 35 

Preedman, W. L., & Madore, B. P. 2010, ARA&A, 48, 
673 

Hanna, C., Mandel, L, & Vousden, W. 2014, ApJ, 784, 

8 

Harry, G. M., & LIGO Scientific Collaboration. 2010, 
Class. Quant. Grav., 27, 084006 

Holz, D. E., & Hughes, S. A. 2005, ApJ, 629, 15 

Hooper, S., Chung, S. K., Luan, J., Blair, D., Chen, Y, 
& Wen, L. 2012, Phys. Rev. D, 86, 024012 


Hu, Y.-M., Raffai, P, Gondan, L., Heng, I. S., 
Kelecsenyi, N., Hendry, M., Marka, Z., & Marka, 
S. 2015, Classical and Quantum Gravity, 32, 
105010 

Kanner, J., Camp, J., Racusin, J., Gehrels, N., (& 
White, D. 2012, ApJ, 759, 22 

Kasen, D., Pernandez, R., & Metzger, B. D. 2015, 
MNRAS, 450, 1777 

Kulkarni, S. R. 2005, ArXiv Astrophysics e-prints, 
astro-ph/0510256 

Li, L.-X., & Paczyhski, B. 1998, ApJ, 507, L59 

Liu, Y, Du, Z., Chung, S. K., Hooper, S., Blair, D., & 
Wen, L. 2012, Classical and Quantum Gravity, 29, 
235018 

Luan, J., Hooper, S., Wen, L., & Chen, Y. 2012, 
Phys. Rev. D, 85, 102002 

MacLeod, C. L., & Hogan, C. J. 2008, Phys. Rev. D, 
77,043512 

Messenger, C., & Read, J. 2012, Physical Review Let¬ 
ters, 108,091101 

Messenger, C., Takami, K., Gossan, S., Rezzolla, L., 
& Sathyaprakash, B. S. 2014, Physical Review X, 
4, 041004 

Metzger, B. D., & Berger, E. 2012, ApJ, 746, 48 

Metzger, B. D., et al. 2010, MNRAS, 406, 2650 

Mohanty, S. D., Marka, S., Rahkola, R., Mukherjee, 
S., Leonor, L, Prey, R., Cannizzo, J., & Camp, J. 
2004, Classical and Quantum Gravity, 21, 765 

Nissanke, S., Holz, D. E., Dalai, N., Hughes, S. A., 
Sievers, J. L., & Hirata, C. M. 2013a, ArXiv e- 
prints, 1307.2638 

Nissanke, S., Holz, D. E., Hughes, S. A., Dalai, N., & 
Sievers, J. L. 2010, ApJ, 725, 496 

Nissanke, S., Kasliwal, M., & Georgieva, A. 2013b, 
ApJ, 767, 124 

Nissanke, S., Sievers, J., Dalai, N., & Holz, D. 2011, 
ApJ, 739, 99 

Nuttall, L. K., & Sutton, P. J. 2010, Phys. Rev. D, 82, 
102002 


8 


Nuttall, L. K., White, D. J., Sutton, P. J., Daw, E. J., 
Dhillon, V. S., Zheng, W., & Akerlof, C. 2013, 
ApJS, 209, 24 

Paczynski, B. 1986, ApJ, 308, L43 

Petiteau, A., Babak, S., & Sesana, A. 2011, ApJ, 732, 
82 

Punturo, M., & et al. 2010, Classical and Quantum 
Gravity, 27, 194002 

Raffai, P, Gondan, L., Heng, I. S., Kelecsenyi, N., 
Logue, J., Marka, Z., & Marka, S. 2013, Classical 
and Quantum Gravity, 30, 155004 

Sathyaprakash, B. S., Schutz, B. E, & Van Den 
Broeck, C. 2010, Classical and Quantum Gravity, 
27, 215006 

Schutz, B. E. 1986, Nature, 323, 310 

—. 2011, Classical and Quantum Gravity, 28, 125023 

Sidery, T., et al. 2014, Phys. Rev. D, 89, 084060 

Singer, L. P, & Price, L. R. 2015, ArXiv e-prints, 
1508.03634 

Singer, L. P, et al. 2014, ApJ, 795, 105 

Soderberg, A. M., et al. 2006, ApJ, 650, 261 

Tanvir, N. R., Levan, A. J., Eruchter, A. S., Hjorth, J., 
Hounsell, R. A., Wiersema, K., & Tunnicliffe, R. L. 
2013, Nature, 500, 547 

Taylor, S. R., & Gair, J. R. 2012, Phys. Rev. D, 86, 
023502 

Taylor, S. R., Gair, J. R., & Mandel, I. 2012, 
Phys. Rev. D, 85, 023535 

The Virgo Collaboration. 2009, technical Report VIR- 
0027A-09 

Wen, L., & Chen, Y. 2010, Phys. Rev. D, 81, 082001 

Wen, L., Pan, X., & Chen, Y. 2008, Journal of Physics 
Conference Series, 122, 012038 

White, D. J., Daw, E. J., & Dhillon, V. S. 2011, Clas¬ 
sical and Quantum Gravity, 28, 085016 

Yang, B., et al. 2015, Nature Communications, 6, 7323 

Zhao, W., van den Broeck, C., Baskaran, D., 8c Li, 
T. G. P. 2011, Phys. Rev. D, 83, 023005 


This 2-column preprint was prepared with the AAS I^TgX macros 


v5.2. 


9 




